Magnetic excitations in vanadium spinels 
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We study magnetic excitations in vanadium spinel oxides AV2O4 (A=Zn, Mg, Cd) using two models: first 
one is a superexchange model for vanadium S = 1 spins, second one includes in addition spin-orbit coupling, 
and crystal anisotropy. We show that the experimentally observed magnetic ordering can be obtained in both 
models, however the orbital ordering is different with and without spin-orbit coupling and crystal anisotropy. 
We demonstrate that this difference strongly affects the spin-wave excitation spectrum above the magnetically 
ordered state, and argue that the neutron measurement of such dispersion is a way to distinguish between the 
two possible orbital orderings in AV2O4. 

PACS numbers: 71.27.+a, 75.30.Mb, 75.20.Hr, 75.10.-b 



I. INTRODUCTION 

Due to geometrical frustration, transition metal spinel ox- 
ides with a general formula AB 2 0^ display a variety of un- 
usual low-temperature properties. The spin dynamics of these 
systems is usually described by a Heisenberg antiferromagnet 
on the pyrochlore lattice. This model is rather peculiar and its 
classical ground state is highly degenerate. Exact degeneracy 
can be lifted by various mechanisms, but the system still pos- 
sesses many competing spin configurations with almost equal 
energies. As a result, when the temperature goes down, the 
system can evolve in a variety of ways: it can remain spin 
liquid down to the lowest temperatures due to quantum fluc- 
tuations, or choose a particular configuration either via the 
strongest order from disorder mechanism or through a struc- 
tural phase transition which lowers the local symmetry of the 
lattice. 

A peculiarity of transition metal spinels is that their mag- 
netic ions often possess also an orbital degree of freedom. 
This extra degree of freedom modulates the spin exchange and 
can at least partially lift the geometrical degeneracy of the un- 
derlying lattice. However, the orbital degrees of freedom in 
geometrically frustrated lattices are by themselves frustrated, 
i.e., many different orbital configurations have the same en- 
ergy. In this situation, the ordering in the orbital sector is cou- 
pled to the ordering in the spin sector, and the selection of 
a true ground state configuration becomes a non-trivial phe- 
nomena. 

Of particular interest is the orbital ordering in the t^g sys- 
tems on the pyrochlore lattice. In this work we study vana- 
dium spinels of the type AV2O4, where A is a divalent ion 
like Cd 2+ , Zn 2+ , or Mg 2+ . In these compounds, magneti- 
cally active V 3+ ions form a pyrochlore lattice and have two 
3d electrons in t2g-orbitals. Due to strong Hund's interac- 
tion, these two electrons form a state with S — 1. All AV2O4 
compounds show qualitatively similar structural and magnetic 
behavior, independently of what the divalent A ion is, and un- 
dergo two phase transitions - a structural one and an antifer- 
romagnetic one. We will be mainly discussing the physics 
of ZnV204. The structural transition occurs at a temperature 
T5 ~ 50 K. 1 Below T5, the lattice shows a tetragonal dis- 



tortion - the vanadium octahedra VC>6 are uniformly flattened 
along the c axis, and the symmetry is lowered from the cubic 
one to IA-i/amd, which is the highest-symmetry tetragonal 
space group for the spinel structure. The antiferromagnetic 
(AFM) transition occurs at a slightly lower temperature Tjv 
of about 40 K. 2,3 This temperature is significantly lower than 
the Curie- Weiss temperature Tew ~ 1000 A', extracted from 
high temperature susceptibility, 4 which underlines the impor- 
tance of geometrical frustration. 

The magnetic structure of AV2O4 spinels at T < Tm was 
first proposed by Niziol 2 and recently confirmed by Reehius 
et al. in Ref. 3. Along the diagonal [1 10]/[110] direction in 
the rry-plane, the ordering is antiferromagnetic H — v while 
along the two other diagonal directions [011]/[011] in the yz- 
plane and [101]/ [101] in the sz-plane, the spin ordering is in 
the form + H v H ... (see Fig.[T). 

At high temperatures T > T$ > Tjv, inelastic neutron scat- 
tering data by Lee et al£ on the powder sample of ZnV2C<4 
still show strong low energy magnetic excitations which form 
a broad peak centered at Q = 1.35 A -1 . This broad peak is 
present also at Tjy < T < Ts, however it becomes asymmet- 
ric and shifts towards a smaller value of Q. The asymmetry of 
the peak further increases in the magnetically ordered phase 
(T < T/v)i£ The value of Q and its temperature evolution 
cannot be explained within a purely spin model, whose fun- 
damental degrees of freedom are antiferromagnetic hexagonal 
spin loops with Q = 1.5 A -1 , which is larger than the exper- 
imental Q. The spin dynamics, however, can be understood if 
one assumes that spin degrees of freedom are affected by the 
orbital degrees of freedom. In particular, the spatial asymme- 
try of the peak in the intermediate phase can be understood 
as the consequence of the fact that at T < Ts, vanadium oc- 
tahedra are flattened and xy-orbital at each site is occupied. 
This leads to a strong antiferromagnetic exchange between 
vanadium spins along xy direction, and, as a result, spin in- 
teractions become effectively one-dimensional. Strong one- 
dimensional spin fluctuations give rise, via spin-orbit cou- 
pling, to fluctuations of the occupations of xz or yz orbitals, 
causing the anisotropy of the neutron peak. 

There were several theoretical attempts to understand the 
nature of the ground state of ZnV^CU. However, although 
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it is widely accepted that orbital degrees of freedom play an 
important role, no consensus is reached yet about the type of 
the orbital ordering (OO) in the ground state. The first at- 
tempt to explain the physics of ZnV204 focused on the spin- 
lattice coupling mechanism, 6 however it didn't explain why 
the structural and spin order occur at different temperatures. 
Tsunetsugu and Motome 7 later addressed this issue and re- 
lated the presence of two separate phase transitions at Tm and 
Ts to the interplay between geometrical frustration and dda 
superexchange (SE) interaction between V-ions. The ground 
state orbital ordering suggested in Ref. |7| consists of stacked 
ab planes with alternating orbital occupations (xy, xz) and 
(xy,yz). Hereafter we label this orbital patterns as ROO. 
They also showed that this ordering of orbitals can partially 
remove magnetic frustration and explain experimentally ob- 
served ordered magnetic structure. 

On the other hand, Tchernyshyov^ pointed out that the 
ground state obtained in Ref. |7| is at odds with x-ray and neu- 
tron diffraction data, because it does not possess the required 
I4i /amd space symmetry. He argued that the spin-orbit cou- 
pling should be included into consideration. He considered 
a purely ionic model in which spin-orbit (SO) coupling plays 
the major role and determines the orbital order in the tetrago- 
nal phase. He suggested the following OO: one electron on 
each site occupies xy orbital, while the second electron is 
spreaded between xz and yz orbitals in such a way as to mini- 
mize the spin-orbit energy at each site. Hereafter we label this 
orbital pattern as COO. 

Recently R. Valenti et al. 9 found in the ab-initio DFT calcu- 
lations that a correct space symmetry can be actually obtained 
within the reasoning of Ref. 7, if one includes into consider- 
ation an additional trigonal distortion. Still, the true ground 
state turns out to be the same as in Ref. la 

The ideas of Refs. and H were combined in the unique 
framework by S. Di Matteo et al. in Ref. 10. They proposed 
to construct a classical ground state phase diagram by consid- 
ering SE interaction and SO coupling on equal footing. They 
demonstrated that the SO coupling is a relevant perturbation 
and favors the states with unquenched orbital momentum for 
any value of the coupling strength. They obtained a variety 
of phases and found that for reasonable values of SE and SO 
couplings, the ground state agrees with the experimentally ob- 
served one. 

In this paper we extend the analysis of Ref. 10 and study 
the low energy excitations in vanadium spinels. We find that 
the magnetic excitation spectrum strongly depends on the type 
of the OO, and that it is qualitatively different for the ground 
states with patterns that consist of real orbitals and those with 
the complex linear combination of orbitals, i.e. complex or- 
bitals. The former ground state is characterized by a quenched 
orbital angular momentum (L = 0), while the later by a un- 
quenched (L 7^ 0) orbital angular momentum. 1 1 The differ- 
ence in the magnetic excitation spectrum arises due to the fact 
that the magnetic moment of the vanadium ion is formed by 
both spin and orbital momentum, and fluctuations of L con- 
tribute to the spectrum of magnetic excitations. We argue that 
the measurement of magnetic excitations in neutron scatter- 
ing experiments can shed light on the nature of the OO in the 



ground state. 

This paper is organized as follows: in Section II we intro- 
duce the model appropriate for the description of the physical 
properties of ZnV204. In Section III we discuss the ground 
state and the magnetic excitations of the system when the or- 
bital angular momentum is quenched. We derive linear spin- 
wave theory of S = 1 moments interacting on the V 3+ py- 
rochlore lattice. In Section IV we discuss ground state and 
magnetic excitations of the system when the orbital angular 
momentum is unquenched. We show that magnetic excita- 
tions for unquenched orbital angular momentum can be de- 
scribed in the framework of the magnetic excitonic model. 
Section V presents the conclusions. Some mathematical de- 
tails are given in the Appendices A and B. 



II. THE MODEL 

The minimal model describing the low energy physics of 
vanadium spinel is given by 

H = H SE + H a + H so . (1) 

The first term describes nearest neighbors (nn) super- 
exchange interactions between vanadium S = 1 spins, arising 
from the virtual excitations djd 2 ; — > d)d^. These interactions 
can be written as: 

H SE = ~ Y«ij> [ J S i ■ Sj + Jl]Oij 

- Y,<ij> J 2 [i - s i • Sj] Oij , 

where i and j are nearest neighbors, Jo = wJ/[l ~ 3?j], 
,h = J[l - v ]/[l ~ 377], J 2 = J[l + ry]/[l + 2ry] are 
coupling constants, J = t 1 j\J\ is the overall energy scale 
(t = "i/Atdda and U\ is the intra-orbital Coulomb repul- 
sion), 77 = Jh /U\ is the normalized Hund's exchange. We 
consider only the largest the hopping term, associated with 
cr-bonding 13 . Such hopping is diagonal and non-zero only 
if the orbitals and the plane in which hopping occurs are of 
the same a(3 type {a(3 = xz, yz, xy). In this case or- 
bital contributions Oij and Oij are expressed in terms of 
projectors P^ a /3 onto the occupied orbital state a(3 at site i 
and j: = P l . a p(l - Pj. a p) + Pj,ap(l - Pi,afi) and 

Oij Pi,a@Pj,a(3- 

To describe the anisotropy and spin-orbit coupling term, we 
use the fact that, when the crystal field splitting between t 2g 
and e g orbitals is large, the t 2g -electrons can be represented 
by an effective orbital angular momentum L' = l. 1 — The 
anisotropy term (the second term of the Hamiltonian (Q])), is 
then given by 

H a = cJ2L7i- (3) 

i 

where c is a constant. This term describes the tetragonal dis- 
tortion in the t^g manifold. We notice, that here for simplicity 
we neglect the trigonal distortion, which is small compared to 
the tetragonal one. 
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The spin-orbit coupling term ( the third term in (fl]i) is given 

by 

H so = -\J2V l -Si, (4) 

i 

where A is SO coupling constant. Note that the true angular 
momentum L is related to an effective one as L ~ —I/. 

The parameters of the model can be estimated from the 
experiments. The spectroscopy data 14 yield the Hund's ex- 
change, Jh — 0.68 eV and Coulomb intra-orbital repulsion 
Ui ~ 6 eV. The estimate of hopping matrix element from x- 
ray photoemission spectroscopy is t ~ —0.35 eV. 1 - 5 , so the 
energy scale is J = t 2 /U\ ~ 20.4 meV. Since the SO cou- 
pling constant is A ~ 13 meV (Ref. 15), we can see that the 
super-exchange and the spin-orbit couplings are comparable 
and, therefore should be treated on equal footing. 



III. REAL ORBITAL ORDER 

A. Ground state 

First we discuss the ground state of the super-exchange 
Hamiltonian alone (Eq. [2}. Without the anisotropy and spin- 
orbit terms the ground state orbital patterns consist of only real 
orbitals: at each site two out of three t 2g orbitals xy, xz, yz, 
are occupied. This type of orbital patterns is called real orbital 
order, as opposed to complex orbital order, when the orbital 
state is formed by a complex superposition of t 2g orbitals in 
such a way, that the gain of spin orbit interaction energy is 
maximized. 

Depending on which orbitals are occupied, one obtains two 
types of interacting bonds and also non-interacting bonds. 
Consider for example a bond in the a/3-plane. If there is 
an electron occupying a(3 orbital only on one site of such 
bond, then the bond, which we label as b\, is weakly ferro- 
magnetic (FM), and is described by H\ 3l = — JoS; ■ Sj — J\. 
If both i and j sites of ij bond are occupied by a(3 elec- 
trons, then the bond, which we label as b 2 , is strongly anti- 
ferromagnetic (AFM). The exchange coupling is then given 
by Hb 2 — — J^il — Si ■ Sj). When neither i nor j site have 
a/3-orbital occupied, the bond is non-interacting. 

One can easily demonstrate 10 that for positive Jo and J2, 
the lowest energy configuration corresponds to the state with 
four ferromagnetic b\ bonds per tetrahedron. There still ex- 
ist two topologically different tetrahedral configurations with 
four b\ bonds, characterized by different OO patterns. One of 
them, the ROO state with OO patterns proposed by Tsunet- 
sugu— , is compatible with experimentally observed magnetic 
structure - it yields AFM chains running in [110] and [110] 
directions (see, FigJTJ- The classical energy per site in the 
ROO state is -Eroo = — 2 J\ — 2 J 2 . However, this state is ac- 
tually not the true ground state of the super-exchange Hamil- 
tonian (0. The other state with 4 b\ bonds (the ROO-I state 
in our notations), in which spins of each tetrahedron form a 
fully collinear up-up-up-down (uuud) state, has the ground 
state energy Eroo-i = —Jo — 2Ji — 2J2, which is lower 




FIG. 1: (Color online) Magnetic ordering consisting of one- 
dimensional antiferromagnetic chains in the xy-plane. Black and 
grey colors correspond to the orbital configuration in the ROO state. 



than that of the ROO state. This is the lowest classical en- 
ergy that one can obtain within the manifold with quenched 
angular momentum at each site 10 . However, the magnetic or- 
dering associated with the ROO-I state is incompatible with 
the experimentally observed one: in ROO-I two neighbor- 
ing spins in xy-plane are ferromagnetically aligned whereas 
the experimentally detected coupling in icy-plane is antiferro- 
magnetic. This discrepancy between the ROO-I ground state 
of the super-exchange Hamiltonian (ffj) and the experimental 
findings demonstrates the necessity of taking into account ad- 
ditional interactions We study the the anisotropy term and the 
spin-orbit coupling in the Section IV. 

One can also consider the model with both the super-exchange 
(SE) and the Jahn-Teller (JT) couplings 17 with a hope that the 
cooperative JT effect, which plays an important role in the 
structural transition, can also stabilize the ROO type of orbital 
order. This model has been studied in Ref. 17 by both the 
mean-field analysis and Monte Carlo simulations. Here we 
assume that the ROO phase can be realized and in the next 
subsection we derive magnetic excitations spectrum for the 
corresponding effective spin model. 



B. Spin waves 

The starting point for the calculation of the magnon excita- 
tion spectrum is the classical Neel ground state with antiferro- 
magnetic spin chains in xy-planes (Fig.Q]). This ground state 
has the magnetic unit cell with eight vanadium spins, which 
we denote as a, b, c, d, d, b, c, d (FigJT}. After averaging the 
orbital operators we can re-write the Hamiltonian (f2]l as 

Hnn — Jxy ^ ' S{ ■ Sj -(- J ^ ' Sj • Sj (5) 

(ij)\\xy {ij)\\xz,yz 

where the first term describes the super-exchange along xy- 
chain with J xy = J 2 , while the second term corresponds to 
the frustrated ferromagnetic inter-chain coupling, J' = Jo- 
To describe the excitation spectrum of such eight-sublattice 
antiferromagnet, we introduce eight boson operators: 
a, b, c, d, a, b, c, d. We employ Holstein-Primakoff transfer- 
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mation, e.g. for up-spins S a and down-spins we have: 

6+6 



S* = S - at a 

s~ 



at a a 



S* b =-S 



(6) 



= aV2S- at a 5,7 = V2S - Wb b 



In the linear 



spin wave approximation we substitute 
V25 in the expressions above. Performing 



V2S- 

Fourier transformation pk = -7= Yli exp _ifcri pi, where N 



is the number of lattice sites belonging to one sublattice, we 
obtain the mean field Hamiltonian for magnons 16 



= (at(k),-a(-k)) 



A(k) 
-B(-k) 



B(k) 
-A(-k) 



a(k) 
at(-k) 



i 7 ) 

where we introduced a(k) = (ak, 6k, Ck, dk, dk, 6k, Ck, dk)- 
The matrices A(k) and B(k) are eight by eight matrices 
whose elements depend on the geometry of the lattice and 
the type of the magnetic ordering, tilde denotes the complex 
conjugation. The explicit expression for matrix elements are 
presented in Appendix A. The quadratic form is diagonalized 
using the generalized Bogoliubov transformation. In the diag- 
onal form, the Hamiltonian takes the form: 

#k,-k = #k : -k + X! A,ik & "k 

n 

+ 5Z A «-k Ck^-k (8) 



where the index n runs from 1 to 8, A„k = A„_k are magnon 
energies and 6„k are linear combinations of boson operators 
belonging to a(k) and a(— k). 

We obtained the spin-wave excitation spectrum by solv- 
ing Eq. ([8]) numerically. The result is presented in Fig. [2] 
along high-symmetry directions of 3D Brillouin zone. We 
used J xy — 18.5 meV, and the ferromagnetic constant J' = 
—0.1 J xy . There are 4 different branches of the spin-wave 
spectrum, each of them is doubly degenerate. Furthermore, 
two of the branches have zero energy over a finite range of 
momenta (see Fig.[2](solid lines)). This so called zero modes 
emerge because the number of ferro- and antiferro-bonds con- 
necting two neighboring antiferromagnetic chains is the same 
(see Fig.[TJ), and the spins forming xy— chains can collectively 
rotate with no change in energy. 

The existence of the zero modes is inconsistent with the 
observation of the magnetic ordering transition at ~ 40K. The 
ordering requires that the zero mode be lifted. The issue is 
what interactions are responsible for the lifting of spin degen- 
eracy. 

A natural first step would be to consider longer range inter- 
actions 17 as these interactions generally remove the degener- 
acy (it happens, e.g., in a Kagome antiferromagnet). We show 
below that it is indeed the case, however, the energy of the 
relevant degeneracy-breaking mode is very small and can not 
explain the magnetic ordering temperature of ~ 40K. 

In Fig[3] we show the interactions which include up to 
third neighbors. There are 6 nearest neighbor interactions 
Jn (Jn = Jxy along xy chain, and J n — J' along xz and 




FIG. 2: (Color online) The spin-wave dispersions obtained along the 
main directions of the 3D BZ. Solid line corresponds to the magnon 
spectrum when J3 = 0, dashed line - to the magnon spectrum with 
the correction due to the third nearest neighbor interaction with J3 = 
0.01 J X y . We have used the following labels for high-symmetry 
points: V = (0, 0, 0), M = (tt/4, tt/4, 0), Y = (0, tt/2, 0), M+ = 
(tt/4, tt/4, tt/4), Y+ = (0, vr/2, tt/4), Z = (0, 0, tt/4). 



yz bonds), 12 second neighbor interactions J nn , and 12 third 
neighbor interactions J nnn . Quite often already inclusion of 
the second neighbor exchange lifts the degeneracy. However, 
here second-neighbor interactions J nn are frustrated and can 
not remove the degeneracy, and, therefore, zero modes in the 
spin-wave spectrum. 17 Thus, one has to include third neighbor 
exchanges. There are two inequivalent sets of third neighbors, 
one obtained by two nearest neighbor steps J nnn and the other 
through the empty space J' nnn (see Fig[3]). 

When only dda hopping is considered, the exchange cou- 
pling through the empty space is zero, J' nnn = 0, and only 
Jnnn interactions contribute. These interactions are antiferro- 
magnetic Jnnn = J3 > 0, and non-zero only if they connect 
sites along the direction corresponding to the symmetry of the 
orbital occupation, (i.e., for orbital occupation a/3, the inter- 
action is nonzero only along a{3 direction). 

The third neighbor interaction is frustrated along xy-chains, 
but it is small compared to nearest neighbor exchange along 
the chain, J3 -C J xy , and can not destroy antiferromagnetic 
ordering along the chain. Along xz and yz directions, J3 are 
not frustrated, and connect parallel antiferromagnetic chains 
located in second neighboring xy planes. The energy scale 
t 4 

for J3 is then -4jg-. This is a very small energy, only about 
one percent of the frustrated ferromagnetic inter-chain cou- 
pling J' ~ -fjf~- This small interaction can not explain the 
magnetic ordering temperature of ~ 40 KJ£ Here we con- 
sider this interaction only qualitatively and assume the value 



of coupling constant J3 = 0.01 J xy . One can see in Fig. [2] 
(dashed lines) that the zero energy modes indeed become dis- 
persive for J3 > 0. 

We would like to note that experimentally it has been 
proven that all magnetic moments are aligned along the in- 
direction. This experimental fact cannot be explained in the 
framework of the SE model, because it is isotropic in a spin 
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FIG. 3: (Color online) J„, J nn , and J„ n n are first-, second- and 
third-neighbor exchange couplings. Second-neighbor interactions 
J nn are frustrated. The third-neighbor exchange coupling through 
the empty space we assume equal to zero, J' nnn = 0; the other cou- 
pling are equal to J n nn = J3, only if they connect sites along the 
direction corresponding to the symmetry of orbital occupation. Red 
and green colors denoting vanadium ions correspond to orbital con- 
figurations in ROO state. 



space. In reality, vanadium spinels likely possess a single ion 
spin anisotropy which aligns spin along the z-direction. The 
anisotropy does affect the magnetic excitation spectrum which 
in its presence acquires a gap. In this case, the zero mode 
will be lifted and the spin frustration will be removed. How- 
ever, since the strength of magnetic anisotropy is not known 
experimentally at the moment, we cannot estimate the magni- 
tude of the anisotropy-induced gap and check whether or not 
this interaction alone can stabilize the ground state magnetic 
structure. We believe that the role of the anisotropy deserves 
further experimental and theoretical investigation. 



IV. COMPLEX ORBITAL ORDER 
A. Ground state 

We now consider the ground state of the system in the pres- 
ence of the SO coupling and the anisotropy term (Eqs. (01 
and (|4|i). For any finite A the SO coupling prefers the orbital 
state with the unquenched effective orbital angular momen- 
tum L' = 1. In such orbital state, one electron at each site 
occupies xy orbital due to the tetragonal distortion, while the 
second electron occupies the complex linear combination of 
xz and yz orbitals. The effective L' should then be parallel to 
the spin magnetic moment in order to minimize the spin-orbit 
energy, i.e. a spin-up site will have L' z = 1 while a spin- 
down site will have L' z = — 1. As we discussed earlier, this 
type of orbital ordering is a COO state suggested first by Tch- 
ernyshyovin Ref. 8. Its energy Ecoo = — l/2[5J2 + 2Ji] — A 
is lower than the energy of the ROO state for a wide range of 
parameters (see the phase diagram in Ref.flol). 

The COO state is characterized by two strong AFM bonds 



per tetrahedra, and its magnetic structure consists of AFM 
chains in xy-planes with the same interaction along the chain 
as in the ROO state. The strengths of the inter-chain cou- 
pling J' are also practically the same for the ROO and the 
COO states 10 , although J' = 1/4[J 2 - 2J ] is antiferromag- 
netic in the COO state, while it is ferromagnetic in the ROO 
state. In this work we assume | J'\/J xy = 0.1 in both orbital 
states. In the COO state the preferred spin direction is fixed by 
anisotropy term (01 via the SO coupling to be along the z-axis, 
and, therefore, spins in xy-chains cannot rotate freely even for 
only nearest-neighbor exchange along the chains. The long 

range order + H ... along diagonal directions in xz and 

yz planes cannot be determined by local interactions, but for 
simplicity we do not include next neighbor hopping terms in 
this part of our calculations. 

B. Magnetic excitons 

We now consider magnetic excitations in the COO state. 
We follow the magnetic exciton model approach of Refs. [T^, 
|2(X which is the extension of the linear spin wave theory for 
systems with unquenched orbital angular moment. 

We consider states with the effective total angular momen- 
tum J = L' + S. Often, energy levels with different J are 
well separated in energy and both J and its z -projection J z 
are good quantum numbers. However, in many transition 
metal oxides the strength of the spin-orbit coupling and the 
super-exchange interaction between localized d-electrons are 
comparable, and atomic energy levels with different values 
of J can cross each other; in such case only J z acts as a good 
quantum number. We show below that magnetic excitations in 
vanadium spinels AV2O4 can be understood as a propagation 
of excitations to states with a given J z through the crystal. 

To proceed, we rewrite the Hamiltonian (Q]) as a sum of a 
single-ion Hamiltonian Hi and the term which describes the 
interaction between two different ions H^- 

H = H 1 +H 2 , (9) 

where 

Hi = Hso + H a + Yli hziS Z i , .jq. 
H-2 = Hse — J2i hziSzi . 

The molecular field part of the exchange interaction acting on 
site i is given by h z — Yl r Z r J r (S z ) r . Z r is the number of r- 
th neighbors, J r is the corresponding exchange constant, and 
(S z ) is the sublattice magnetization. 

First, we diagonalize the single ion Hamiltonian, H\, in 
the molecular field approximation. It is convenient to express 
the eigenfunction \ J Z ) for the states splitted by the spin-orbit 
interaction as linear combinations of the unperturbed eigen- 
functions of L' z and S z . Then H\ can be represented as a 
block 9x9 matrix in the subspace of \L' Z , S z ) as follows: 

H X {J Z = ±2)=c-\±h z (11) 
F l( J,=±l)=(_ C A (12) 
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5 10 15 20 25 30 35 40 45 50 
h z [meV] 

FIG. 4: (Color online) Energy levels of V 3+ ions in AV2O4 as func- 
tions of molecular field h z . J z — —2 is a ground state for all values 
of h z (solid blue bold line). The transitions are possible only to ex- 
cited states with J z = — 1 (solid red thin lines). 



fli(J, = 0) 




(13) 



where J z — 2, J z — —2, J z — 1, J z = — 1 and J z = 
are represented in the basis |1, 1), | - 1, -1), (|1,0), |0, 1}), 
(| - 1,0), |0,-1)), and (|1,-1), |0,0), | - 1,1)), respec- 
tively. Diagonalizing these matrices, we obtain eigenvalues 
and eigenvectors of Hi as functions of the molecular field h z . 

When the tetrahedra are flattened as observed in the exper- 
iments, i.e. when c < 0, for any strength of the molecular 
field h z the ground state of Hi has J z — ±2 (the sign of 
J z depends on the sign of h z ). For definiteness, we consider 
spin-down sublattice (h z > 0). The resulting energy levels are 
presented in the Figj4] The magnetic excitations of the local 
Hamiltonian Hi are transitions from the ground state to eight 
excited states of a single ion. These excitations can be de- 
scribed by boson operators p\ - each pseudoboson describes 
the transition from the ground state |0) to the excited state \ v). 
The single ion Hamiltonian Hi is diagonal in terms of pseu- 
doboson operators: 



(14) 



i V — \ 



where e v is the energy difference between the excited state \ v) 
and the ground state of Hi : 



E v — E . 



(15) 



We next introduce the representation for spin operators in 
terms of pseudobosons p\. The representation for spin S p on 
the sublattice p can be written as: 

s* = (o\s,\o) + yL 1 (v\s,\oM+pv)+ 16 , 
El=i((v\s*\v)-(o\s z \o)Hp v 



S p 



(17) 



If the ground state of Hi corresponds to J z = — 2, it follows 
from ( fT71 ) that only states with J z = — 1 contribute to Sp. 
There are two such states, |1) and |2). In the basis of \L' Z , S z ) 
they have the following structure: 



|i) = VT~^\ - 1,0) +o|o,-i) 

- 2 |0, 1> 



|2) = a|-l,0)-Vl 



or 



(18) 



where we denote a as the weight of the state with L' z = 0, 
and S z = —1 in the |1) state. The partial weight of orbital 
and spin contributions is determined by the competition be- 
tween the spin-orbit coupling, the anisotropy energy and the 
molecular exchange field. In Fig. [5] we present a dependence 
of a on the molecular exchange field h z at different values 
of the tetragonal field c, keeping the SO coupling constant A 
fixed. For all values of tetragonal field parameter a shows 
similar field dependence: it increases with the increase of h z . 
This happens because the molecular field h z , acting on the 
spin Si, is effectively created only by the nearest neighbors 
on the xy-chain, which are antiferromagnetically aligned to 
Si. Therefore, when h z increases, the flipping of the spin 
costs more energy and the transitions with AS Z — 1 are sup- 
pressed. The low lying excitations become more of an orbital 
character, which can be seen in the increase of the weight of 
the transitions with AL' Z = 1. 

As only two excited states of an ion are relevant to the spin- 
wave analysis, we consider the local excitations described by 
two pseudoboson operators p\ and p\, which take an ion from 
the ground state to states |1) and |2), respectively. The explicit 
form of extended Holstein-Primakoff transformation for spin 
S p in terms p\ and P2 is given by: 



S*= (0|£,|0)+ 

«1|5,|1) - {Q\S z \Q))pW 
{{2\S Z \2) - <0|S,|0))i4p2 



5± = < llS^lO > p\+ < 2|5 ± |0 > 

< 015*11 > P i+ < 01^12 >p 2 



pi- 



rn 



(20) 



It is also useful to rewrite these expressions, ( fT9] l and d20b . 
using the definition of a (for spin-down): 



-1 + (1 - a 2 )p\pi 



P2P2 



Sp 



y/(l - a 2 )2p\ + u\[2p\ 



S~ = ^/{l-a 2 )2pi+aV2p 2 



(21) 
(22) 
(23) 



The interactions between localized excitations are de- 
scribed by i?2- The excitation spectrum is obtained in a simi- 
lar way as in the spin wave analysis for the ROO state. After 
diagonalization, the total Hamiltonian H = Hi + H2 can be 
written in the same form as Eqs.dSJ, however now the index 
n runs from 1 to 16 and the modes have complex spin-orbital 
character. The details of calculation are given in Appendix B. 

The numerically calculated magnetic excitation spectrum 
for the COO state is presented in Fig|6] For comparison, 
we also plotted there the spin-wave dispersion for the ROO 
state. We used the following parameters: J xy — 18.5 meV, 
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FIG. 5: (Color online) a as a function of molecular field h z for dif- 
ferent values of the tetragonal field c. 




FIG. 6: Magnetic excitation spectrum. Red solid lines correspond 
to magnetic excitations in case of unquenched orbital angular mo- 
mentum, 1^0, obtained in the framework of the magnetic exciton 
model (MEM). We use the following parameters: J xy = 18.5 meV, 
J' = 0.1J X y, c = —1 meV, and A = 13 meV. The black dashed 
lines correspond to pure spin waves, L = 0, obtained in the linear 
spin wave approximation (LSWA).The spectra are plotted along the 
same high symmetry directions as in Fig[2] 



J' = 0.1 J™, c = — 1 meV, and A = 13 meV (third neighbors 
interactions are not included here). As we see from Fig|6] the 
excitation spectrum consists of 8 different branches, each of 
them is doubly degenerate. The excitations may be divided 
into two groups: four low lying branches with rather small 
gap A Q ~ 3 meV, and four optical branches, with the gap 
A Q ~ 35 meV. As A a -C A a , the lying branches are quasi- 
acoustic modes. As we see in Figj6j the lowest mode along 
the direction T = (0,0,0) -> M = (tt/4, tt/4, 0) is the 
dispersionless mode. This flat mode is a lifted "zero energy 
mode" of the ROO state. Here the zero mode is lifted due to 
the combined effect of the anisotropy term and the SO cou- 
pling. The tetragonal distortion favors the non-zero value of 
the ^-component of orbital angular momentum, which in turn, 
selects the local spin quantization axis in such a way, that the 
z component of spin is non-zero. We found that the gap, A Q 
is mainly determined by the anisotropy term and only weakly 
depends on the SO coupling (see FigO. 

However, we note once again, that, in principle, the lifting 




01 [meV] 

FIG. 7: (Color online) Anisotropy A a (dashed line) and optical A 
(solid line) gaps as a function of spin-orbit coupling constant A. 



of zero energy mode can appear also in ROO state if single ion 
magnetic anisotropy is taken into account. We caution that the 
size of these two different anisotropy gaps can be of the same 
order, and therefore, it will be rather difficult to distinguish 
between them. 

As a remark, we would like to mention, that such lifted zero 
energy modes can be detected by an inelastic neutron scatter- 
ing, and have been indeed recently observed in the frustrated 
Kagome lattice antiferromagnet KFe3(OH)6(S04)2^- 

There is another significant difference between the mag- 
netic excitation spectrum for COO and ROO states: as one 
can see from Fig|6j the bandwidth of quasi-acoustic modes in 
the COO state is strongly reduced compared to the bandwidth 
of pure spin waves in the ROO state. This reduction is the ef- 
fect of the mixing between orbital and pure spin excitations in 
the spectra of the COO state. Pure orbital excitations are non- 
dispersive, because they come form local interactions. The 
presence of the orbital component in the spin-wave spectrum 
then obviously leads to the reduction of the bandwidth. 

We believe that this large reduction of the bandwidth can be 
seen experimentally, in the neutron scattering measurements 
of the magnetic excitation spectrum. Such experimental re- 
sults would discriminate between two types of orbital order- 
ing, COO and ROO. Unfortunately, at present these exper- 
iments are difficult to perform because to measure the full 
spin-wave spectrum one needs single crystals, which are still 
not available, as far as we know. 

Let us now discuss the optical branches. These modes ex- 
ist only in the COO state, and they arise from the hybridiza- 
tion between orbital angular momentum non-dispersive levels 
and dispersive spin branches due to the spin-orbit coupling. 
In FigO we show how the optical gap, A Q , depends on the 
strength of the spin-orbit coupling. One can see that at A = 0, 
A Q is zero, but it rapidly increases with increasing A. 

To summarize, the magnetic excitation spectrum for the 
COO state has two separated branches: quasi acoustic modes 
and optical modes. Both manifolds are gapped, but the gaps, 
A a and A D have different origins: A a is set by the anisotropy 
term, while A Q is set by the spin-orbit coupling and is much 
larger. 

These two gaps could be, in principle, determined by in- 
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elastic neutron or Raman scattering even in powder samples 
of ZnV 2 C>4. They are also easily distinguishable from the 
well-known Haldane gap, Ah, which is a characteristic fea- 
ture of the antiferromagnetic S = 1 chains. Ah, A a and A Q 
have different temperature dependence: the optical gap, A G , 
does not depend on the temperature, as it is determined by the 
relativistic spin-orbit interaction; the anisotropy gap, A a , is 
non-zero only below the temperature of the structural transi- 
tion, T < T s ; and the Haldane gap should disappear below 
Tn, when a long range antiferromagnetic ordering emerges. 

At all temperatures, the optical gap, A Q , is the largest, so 
we compare only the anisotropy and the Haldane gaps. At the 
lowest temperatures, T < Tjy, the Haldane gap vanishes, and 
the lowest mode will have a gap equal to an anisotropy gap 
A a . At intermediate temperatures, Tjv < T < Ts, the spin 
excitations are also gapped. In this temperature range, the gap 
is the sum of the anisotropy gap, A a , and the Haldane gap, 
Ah- The magnitude of the well-developed Haldane gap is of 
the order of 0.4 J xy , and is compared to Ts, hence it is very 
likely that the gap at T N < T < T s will be larger then at 
low temperatures. This behavior is exotic and, as far as we 
know, has not been yet observed in S = 1 spin-chain systems. 
At high temperatures, T > Ts, the sharp gap is washed out 
by thermal fluctuations, but the spectrum can still divided into 
acoustic and optical branches. 



V. CONCLUSION 

We presented in this paper a detail analysis of the mag- 
netic excitations of vanadium spinels, whose low-temperature 
tetragonal phase can be modeled to a high accuracy by one- 
dimensional spin chains with weak inter-chain interaction. 
The formation of antiferromagnetic spin chains on the highly 
frustrated pyrochlore lattice is by itself non-trivial phenom- 
ena. This can happen only because vanadium ions also pos- 
sess an orbital degree of freedom, and the orbital modulation 
of the spin exchange partially lifts the geometrical degeneracy 
of the underlying lattice. 

We considered two different ground states: i) the one with 
the real orbital ordering, ROO, and ii) the one with the com- 
plex orbital ordering, COO. We found that the excitation spec- 
tra in these two cases are qualitatively different. The spec- 
trum for the COO state consists of low lying quasi-acoustic 
modes with small anisotropy gap and optical branches with 
zone-center gap determined by the spin-orbit coupling. The 
spectrum for the ROO state has only quasi acoustic modes. 
Within the superesxchange model considered in the present 
study the spectrum is gapless, however, in reality we expect 
an anisotropy gap also for this state. 

The bandwidth of the quasi-acoustic modes in the COO 
state is strongly reduced compared to the ones for the ROO 
state, due to contributions from orbital L-modes. Because 
the spectra are so different, we argue that an effective way 
to determine experimentally the symmetry of orbital ordering 
in vanadium spinels is to measure their magnetic excitation 
spectrum. 

We gratefully acknowledge discussions with R Fulde, A. 
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A. Appendix A 

In this appendix we present the expressions for the matrices 
A(k) and B(k) in the Hamiltonian ||7). 

The diagonal elements of the matrix A(k) are given by 

App(k) = 2J xy + 2J 3 cos2(fe E - k y ), p= 1,3,5,7 
App(k) = 2J xy + 2J 3 cos2(k x + k y ) 1 p = 2,4,6,8 

Nonzero matrix elements of A(k) are given by 



A M (k) 


= J' 


e 


— i(kx— k z ) 


Ai 6 (k) 


= J' 


e 


~ ^(. k y ~\~k z ) 


A 23 (k) 


= J' 


e 


-i(k x + k z ) 


A 25 (k) 


= J' 


e 


{k v +k z ) 


A 38 (k) 


= J' 


e 


(ky-k z ) 


A 47 (k) 


= J' 


e 




A 58 (k) 


= J' 


e 


-i(kj;—kz) 


A 67 (k) 


= J' 


e 


-i(k x + k z ) 


Apg(k) 


= A 


qp 


(k) 



Nonzero matrix elements of B(k) are given by 



B13 


£k) 


— 2J X yCOs(k x 


ky) 


B24 




— 2J X y cos(k x + 


ky) 


B57 




= 2J X y cos(k x 


ky) 


B68 


,k) 


= 2J X yCOs(k x + 


ky) 


B12 




= J 1 e i( - k v+ k ') 




Bis 




= j'e^C**-**) 




B27 




= J' e i(k x +k z ) 




B34 


;k) 


_ j' e -i(.ky-kg) 




B36 


;k) 


= J' e - ' i ( fe *+ fe *) 




B45 




= J' e ~ i ( k *~ k *) 




B56 




= J' e i( - ky+kz ^ 




B78 




— j' e -i(k y - k *) 




B15 




= 2 J3 cos 2(k x - 


k z ) 


B26 




= 2 J3 cos 2(k y — 


k z ) 


B37 


;k) 


= 2 J3 cos 2(k x -+ 


k z ) 


B48 




— 2J 3 cos2(fc y - 


k z ) 
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B P9 (k) = B gp (k) 



where 



B. Appendix B 



A pq (k) 



A p9 (k) p ^ q 
p = q 



(25) 



We next consider the modifications of the matrices for the 
exciton model. 

In the exciton model, the base object is an enlarged set of 
pseudoboson operators a(k) = (ai(k), a2(k)), whose com- 
ponents ai (k) and &2 (k) describe transitions to first and sec- 
ond excited levels. These two vectors are analogous to a(k) 
for pure spin wave model. As a result, the matrix A(k) is en- 
larged and becomes a 16 x 16 matrix A(k). Its components 
are given by 



A(k) 



(1 - a 2 )A(k) 
aVT- 



A(k) 



ayl 
a 2 A{k) 



a 2 A{k) 

£2 



(24) 



Here p and q runs from 1 to 8. The diagonal matrix elements 
of A pq (k) are zero (we consider only the case J3 = 0), be- 
cause the diagonal contribution is akeady included in e v . 

Similarly, we can obtain the expression for the matrix 
B(k): 



(1- 



a 2 )B(k) aVl - a 2 B(k) 



* 2 B(k) 



<B(k) 



(26) 
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